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Abstract 

We use the method of color effective Hamiltonian to study the properties of states in which a gluonic subsystem forms 
a color singlet, and we will study the possibility that such a subsystem hadronizes as a separate unit. A parton system 
can normally be subdivided into singlet subsystems in many different ways, and one problem arises from the fact that 
the corresponding states are not orthogonal. We show that if only contributions of order l/N^ are included, the 

. problem is greatly simplified. Only a very limited number of states are possible, and we present an orthogonalization 
procedure for these states. The result is simple and intuitive and could give an estimate of the possibility to produce 

. color separated gluonic subsystems, if no dynamical effects are important. We also study with a simple MC the 
possibility that configurations which correspond to "short strings" are dynamically favored. The advantage of our 
C$ \ approach over more elaborate models is its simplicity, which makes it easier to estimate color reconnection effects in 

■ reactions which are more complicated than the relatively simple e + e~ annihilation. 
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I. INTRODUCTION 



Hadronic processes in high energy reactions are generally formulated in terms of two separate phases, one initial 
perturbative phase described by quarks and gluons followed by a nonpertubative hadronization phase. In the first phase 
the probability to produce a specific parton state, with specified colors, can be calculated by perturbative quantum 
chromodynamics(PQCD). The transition from a partonic state to a hadronic one is a nonpertubative phenomenon 
and can only be described by hadronization models. Hence there is an intermediate interface between the perturbative 
stable partonic state and the hadronic one. To describe this interface it is not enough to know the momenta and colors 
of all partons, we also have to know how they are connected. Each color charge must be associated with a partner 
O | anticharge. Therefore the interface is normally described by a string or a cluster chain. This intermediate interface is 
needed not only to determine the hadronization, it is also needed to specify the termination of the perturbative cascade. 
In the HERWIG [QJ|] and ARIADNE |j] cascades the termination is specified by a local transverse momentum and 
in the PYTHIA/JETSET H) cascade by the virtual mass of the parent. In all cases an ordering has to be known. We 
want to stress that in a perturbative calculation, which includes interference effects, the necessary parent-daughter 
relation cannot be specified from perturbative QCD. Instead it depends on the nonperturbative, soft, properties of the 
strong interaction. This implies that also if one assumes e.g. local parton hadron duality (LPHD), where the number 
of hadrons is directly related to the number of partons, a calculation of the multiplicity must rely on assumptions 
about the cutoff which cannot be specified from perturbative QCD alone. (These features have been more thoroughly 
discussed in Ref. ||.) 

This problem does not appear in the large N c limit. When N c — > oo there is a single color ordering, which is the 
only possible one. Planar Feynman diagrams dominate and many interference terms vanish. In this limit we also find 
from perturbative QCD that generally partons are combined so as to produce a rather "short" string or cluster chain. 
The process e + e~ — > qq + ng corresponds to a specific chain stretched by the gluons. In the perturbative phase, the 
only possibility to split the system into two or more separate subsystems, is by the process g — > qq. Thus all strings 
or chains have a quark or an antiquark in the ends, and it is not possible to isolate a purely gluonic subsystem. 

For finite N c , there can be many different ways to connect the color charges. This ambiguity is connected with the 
possibility for two gluons to have identical colors, and is therefore generally of order l/N^. New types of states are 
possible compared to the large N c limit. Thus for the e + e~ annihilation process it is possible to have a purely gluonic 
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subsystem, which can form a color singlet and therefore hadronize as separated from the remainder. (This possibility 
was to our knowledge first discussed in Ref. |(|.) 

We are then facing the question: what decides the hadronization in a situation when several different color connec- 
tions are possible: 

• Does Nature make a random choice among the different possibilities? This means that the problem is determined 
purely by the SU(3) group structure. 

• Are some states dynamically favoured? If this is the case the probability for a specific configuration should 
depend on the combination of a group factor and a dynamic factor. 

• If dynamics is important, are the properties in momentum space more essential, (e. g. favouring states which 
correspond to "short strings" Ref. 0) or are the properties in coordinate space more important (as proposed e. 
g. in Ref. |||)? 

We can isolate three different problems: 

f . The SU(3) group structure - How many possibilities are there, and what are their relative weights? 

2. What are the dynamical properties of the alternative states allowed by the group structure? 

3. What observables can be used to solve the problems? 

In the standard models ||l|,^,^T) for parton cascades and hadronization, the only effect of finite N c originates 
from the coupling of a gluon to quarks, which is given by Cp = hN c (l — 1/-V?) and thus reduced by a factor 
(1 — I/-V?) compared to the triple gluon coupling. Only contributions from planar Feynman diagrams are included 
and a set of color suppressed effects are neglected. This implies that in these models, just as in the large N c limit, an 
e + e~ — > qq + ng event can split into two separate strings or chains only via the process g — > qq. Thus all strings or 
chains have a quark or an antiquark in the ends, and it is not possible to isolate a purely gluonic subsystem. 

As mentioned above, the state intermediate between parton states and hadron states can be described by a string 
or cluster chain configuration. If this state contains a closed loop formed by gluons, we call it a color separate(CS) 
gluonic sub-singlet. 

It is clear that to estimate the production rate of CS states, it is not enough to calculate the probability to find a 
gluonic subset which forms color singlet. For a state with many gluons this probability would always be one, and a 
gluonic singlet subsystem could also be formed in many different ways. 

Due to their non-perturbative nature, there is no way so far to determine the probability of string states with CS 
gluonic sub-singlets. A set of models have been presented, which include nonconvcntional color configurations, where 
the partons can be "reconnected" compared to the "standard" configuration [[^]l2]-|i"4|] . Such a model may easily 
appear as a black box in the sense that it is difficult to judge what property of the model is responsible for a particular 
observable effect. This paper is an attempt to make it easier to understand the features of CS states with gluonic 
sub-singlets. Here we use the method of the color effective Hamiltonian which is constructed from PQCD invariant 
amplitude p"5| to study CS states corresponding to a particular string configuration from our available knowledge of 
PQCD. 

In section II- IV, we study the effects of the group structure using this method. We try to estimate the probability 
to produce a color singlet with gluonic sub-singlets(CSGS), disregarding any dynamic effects which might disfavor its 
hadronization as a separate unit. For the parton states produced in PQCD the production amplitude can be written 
as a sum of terms, where each term contains a trace in color space, which corresponds to a definite color ordering 
of the gluons. We call these parton states \hi). We also want to study a different set of states, called singlet chain 
states, for which the color of a gluon (or quark) and the anticolor of a neighbouring gluon (or antiquark) form a color 
singlet. These states can form a single chain (the states \ fi)) or they can form two or more separate chains with at 
least one isolated gluonic loop (the states \si)). In the large N c limit the states \hi) are identical to the states 
and they correspond directly to string states with the same color ordering. For finite N c this is no longer the case. 
It is not possible to directly associate the singlet chain states, \fi) and \si), to definite string configurations. One 
reason is that the singlet chain states are not orthogonal. They can, however, as will be discussed in section IV, be 
adjusted into an orthogonal set by corrections of order l/N^. We also show that to leading order in l/N% the number 
of possible states \si) is significantly reduced, which greatly simplifies the result. The states obtained in this way may 
be associated with specified string (or chain) configurations. Here we propose a first estimate for the group weight 
corresponding to a specific string configuration obtained by projecting these string candidates upon a specific \hi). 

In section V we present some results from a very simple " reconnection" model based on the rules about leading 
CSGSs which we obtain in the earlier section. The comparison between global properties of events with and without 
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CS states are given. By comparing with more elaborate models for e + e~ annihilation one may hope to clarify if some 
simple regularities are most important for the result. The simple model might then be used to estimate reconnection 
effects in more complicated reactions. 

II. COLOR EFFECTIVE HAMILTONIAN 

The QCD Lagrangian is: 



where A™ (a = 1, • • ■ , 8) are eight gluon fields, q is the quark field with three colors (R, Y, B); field tensor F" is given 
by 

F% v = d^AZ - d v Al + gf a ^A^AZ 

and T a — A Q /2 where A Q is the Gell-Mann matrix for SU(3). We can write the QCD Lagrangian in another form by 
combining Gell-Mann matrix to a set of new ones: 

{T12 = ^(Ai +«A 2 ),r 2 i = 7-^,713 = i(A 4 +iA 5 ),r 3 i = t^,t 23 = ^(A 6 + iA 7 ),T 32 = r£ 3 , A 3 , A 8 } 

Matrix nj raises color j to color i, and its Hermite conjugate = has the opposite function. Accordingly we 
define eight new gluon fields 

x l = ^4 > K = ^4 

Xf = ^{Al + iAl),Xf = l,{Al-iAl) 
Xf = ^{At+iAl),Xf = l={Al-iAl) 

X? = ±(At + iAl),X? = ±(4-iAl) 

Hence the quark-gluon coupling part of the QCD Lagrangian can be written as: 

C int = gq^A«T a q 

= -fetfV?Xl + \*X*)q (1) 
+ ^[qr(T2iXl 1 + T^Xf + n 2 Xf)q + h.c] 

where X^ couples to color RR and YY ; X® to RR ,YY and BB; X 21 to RY ; X 12 to YR; and so on. So we write 
the QCD Lagrangian in a form with the quark-gluon interaction term showing clear color significance. This triggers 
our following color treatment for multigluon processes. 

We can construct from PQCD a strict formulation to calculate the cross section of color singlets of a multiparton 



system at the tree level 15 . For the process e + e~ — > qq + ng, the essential part of the formulation is to exploit the 
color effective Hamiltonian H c to compute the amplitude (/|iJ c |0) of a certain color state |/). The color effective 
Hamiltonian H c is found from the invariant amplitude MS 



P 



where a/b is the color index for quark/ antiquark; a%, a 2 , • • • , cc n are those for n gluons with a u = 1, • • • , 8 (for u — 
l,---,n); the summation is over all permutations of (1, 2, • • • ,n); D p = D(q,q,gpn\,gp/2^,- ■ ■ ,gpt n )) is the function 
of parton momenta where momentum indices are suppressed, and P denotes a certain permutation of (1,2, • • ■ ,n); 
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(Z>a P(1)T a Pl2) . . .T Qp (")) Qb is the a-th row and the 6-th column element of the matrix (T ap wT ap m ■ ■ ■ T ,Qp ("> ) . The 
form of the invariant amplitude can be derived from the ST/ (3) color structure of QCD. A modern and convenient 
treatment was developed by Berends and Giele |l6[| . They proposed a recursive method from which eq.(|J) can be 
readily obtained after including all Feynman diagrams at the tree level. Then H c is built from eq.(||) as follows: 

= Ep(l/V2)»Tr(Qt G t 3(i)G t 3(2) ... G t 3(n)) ^P " ' (3) 

where the repetition of two subscripts represents summation(we use this convention unless explicitly specified); ty' a = 

(itf , yt, B^) is the color creation operator for quark; "f^ = (R},Y ,B)is the anticolor creation operator for antiquark; 
(Q )a = * 1S the nonet tensor operator; For gluon u, G' u , the gluon's color octet operator, is defined by: 

Gl = (l/V2)A*^"t = ^f ua _ l ES{aj b) ^ ua = G 'i -\e-1 u (4) 

where GJf = #t f *L 1« = <Ka, and *L are C^j^)^)« and (R\Y\B^) U respectively; E is the unit 

3x3 matrix; 5(a, b) is Kronecker symbol: 5(a, b) = { J' ^ _^ ^ ; A""''' is the a u -th gluon color operator whose meaning 

resembles that of field. Note that we use the same symbol A to denote these color operators as the gluon field in 
the former paragraph and this does not cause ambiguity because we no longer use the gluon field in the rest of the 
paper. 

The color effective Hamiltonian is another expression of the S matrix, so it is not Hermitian. The validity of H c 
can be verified by the calculation of the matrix element for the process e + e~ — ► qq + ng. For the color initial state 

|0) and the final state |/) = ty a '^ b A^A^ 2 ■ ■ ■ A""^, summing over the color indices a'/b' for the quark/antiquark 

and those a'^, a' 2 , • • • , ot' n for n gluons, and carrying out all phase space integrals, we obtain: 

/ dfl Y, s K/l He |0) | 2 = JdQ (0| H\H C |0) 

= JdQ M a Q b ia2 "' Q " • (M a Q b ia2 "' Q ")* (5) 
= <T tr ee(e + e~ — > qq + ng) 

where J dQ, denotes phase space integrals including all kinematic factors. We can see that the calculation of the total 
cross section via H c returns to the ordinary form. 

A color state of the multiparton system qq + ng is composed of color charges of q, q and n gluons. It belongs to the 
color space: 

3^ 3| <g) 8 X 8 2 <g) • • • 8„ (6) 

or a little larger space: 

3 9 3^ 3 X 3* 3 2 5g • • • 3„ 3* (7) 

The effects of unphysical " singlet-gluon" state brought by enlarging the color space 8„ to 3 n 3* (= l n 8„) in (Q) 
can be eliminated in calculating the matrix element of H C1 i.e. projecting the state to H c |0), because 

(f u ■ the rest of the color state\ H c |0) = 

where u denotes a gluon. In other words, (/ | H c | 0) = if |/) contains the singlet composed of the color and its 
anticolor of a single gluon. There are many ways of reducing the color space (|^) or (Q) . If a complete and orthogonal 
singlet set is denoted by {\fk), k = 1, 2, • • •}, we have: 

|A) (A- 1 = 1 (completeness), (A I A') = hk> [orthogonality) (8) 

and 

E fc = / dn £ fc I <AI He |o) | 2 = / dtt <o| Ft | /fe) (/fe | # c | ) 

= /dfi M™""™ • (ilC """**)* = ^ee(e+e- - ?5 + n fl ) 1 j 

This property implies that the sum of the cross sections over all color singlets in a complete and orthogonal singlet 
set for the system (qq + ng) is equal to the total cross section at the tree level. 
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III. EXAMPLES OF COMPLETE AND ORTHOGONAL SINGLET SETS 



A. e + e — > qqg 

For the process with only one gluon in the final state, i.e. e + e~ — ► qqg\, the color space is: 

3^0 3^0 8! (10) 

The enlarged color space is: 

3 g <8> 3± ® 3j 3J (11) 

There is only one way of reducing the space (|Io|). It corresponds to a unique singlet of the system qqg. This singlet 
is \Tr(QGi)), where {Q) a is the color nonet tensor of the quark pair and can be written as the irreducible octet tensor 
(Q')a plus its trace 1q = Tr(Q): Q = Q' + \1qE. Thus the complete singlet set is composed of only one singlet: 

{\Tr{QG 1 ))} = {\Tr{Q'G 1 ))} (12) 

where the equality is due to the fact that the gluon octet tensor is traceless, i.e. Tr{G\) = 0. 

There are many ways of reducing the space ([ll]). Here we choose the following one that can lead to the singlet 
chain state: 

(3, 3D (3± 3 X ) = (l 9l © 8 ql ) (Ii, © 8 1? ) (13) 
Corresponding to this reduction is the complete and orthogonal singlet set as follows: 

{i|l,il 1? ),i|8 9 i®8i ? )} (14) 

where l uv = 8{a u , by^a^ and u and u denote the parton, a u /b v is the color /anticolor index for parton u/u; and 

18™ ® 8 st ) = Tr(Gl v Gt t ) |0) = |Tr(G w G rt )) 

where u, v, s, f are parton labels and {Gl v ) b a = **t*t - |M(a, 6)*^*^. The first singlet of @ is defined as the 
singlet chain state. Generally a singlet chain state starts from the quark end, connects n gluons one by one in an 
order and ends at the antiquark end. Each piece of the chain state is a color singlet formed by the color of one parton 
and the anticolor of the next parton in this order. The second singlet of ( |i"4| ) is the orthogonal one of the singlet chain 
state and is obtained from the reduction 8 ql eg) 8-^. 



The color space is: 



The enlarged color space is: 



B. e+e qqgg 



3 9 ® 3£ ® <g> 8 2 (15) 



3 g ® 3^ <g> 3j ® 3J ® 3 2 ® 3; (16) 
Let us look at one reduction way for space (|l5|): 

(3, 3±) (§! ® 8 2 ) = (l g5 8 8^) ® & 8 2 ) (17) 
The complete and orthogonal singlet set corresponding to this reduction is: 

{(2V6)- 1 |l^(8i ® 82)) , (3/80) 1 / 2 |8, 5 {81 , 8 2 }) , (48)- 1 / 2 |8 ?5 [8 X , 8 2 ])} (18) 
where the symmetric octet state |{8i ,8 2 )} and the anti-symmetric one |[8i ,8 2 ]) are defined as: 

2 



|{8i,8 2 }) 



8{a u b 2 ){Gi) b ai (G a )* + 8(a u 6 2 )(G 2 )g 2 (G 1 )« 1 - -*(o, b)Tr(GtG 2 ) 
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|[8 1 ,8 2 ]) = |J(a 1 ,6 2 )(G 1 )^ i (G 2 )^-«5(a 2) 6 1 )(G 2 )^(C? 1 )^) 
Another example is the reduction of the enlarged space ([!(]) which can give the singlet chain state: 

= © 8,1 ) ® (li2 © 8ia) ® (1 2? © 8 2 ^) 1 y) 

The complete and orthogonal singlet set corresponding to this reduction is: 

{(3V3)- 1 |l gl l la l33> , (2V6)- 1 |l,i(8i2 ® 8 2 q)) , (2 v / 6)- 1 |li 2 (8,i (8 8 2? )> , . . 

(Sy^)" 1 |l2?(8 9 i ® 8i 2 )> , (3/80) 1 / 2 |8 25 ® {8 g i , 8i 2 }> , (48)^/2 |g 2 _ g, [ 8gl , 8l2 ] )} I J 

In (|2C|), the first singlet is the singlet chain state which has the parton configuration: (qgig^q) ■ One can interchange 
gluon labels 1 and 2 in ( [l9| ) to obtain another singlet set which includes the configuration (qg 2 giq). 

IV. CSGS IN QQ + NG SYSTEM 

As we see in the second section, the color effective Hamiltonian H c contains all perturbative dynamics relevant to 
color interaction. H c is actually the color S matrix. We may define the \in) state as we do in quantum scattering 
theory: 

\in) = \H e ) = H c |0) (21) 

We call the state \H C ) the physical color state. We can calculate the color matrix element of any color state |/) by 
projecting the state to the physical one: (/ | H c ). According to (||), the physical state is written as: 

\H C ) - ^(1/V2)"£ P \Tr(QG P{1) G P{2) ■ ■ ■ G P[n) )) (22) 
p 

where there are n\ trace terms. Each of them is associated with one order of gluons: P(l,2,- • ■ ,n — l,n), and it 
corresponds to a string connection of the same order at the large N c limit because they are orthogonal to each other. 
But with iV c = 3, different trace terms overlap and no exact correspondence exists between a trace term with a 
specific gluon order in the physical color state at the partonic level and the string connection of the same order at 
the hadronic one. For each trace term, the color configuration is not uniquely determined and there are ingredients 
of color singlet chain states and CSGSs. In this section, we will discuss some properties of the CSGS and show how 
to estimate their probability from PQCD. 

As the first example, let us look at the simplest case: e + e~ — * qq\g\g 2 - The physical color state for this process is: 

\H C ) = X - \Tr(QGiG 2 )) D 12 + X - \Tr{QG 2 G l )) D 21 (23) 

According JIs|), the only CSGS for the system qqg\g 2 is: 



l s i) = ^7^=ll T KGiG 2 )l^) 

=(iii 2 i 2 ii™)-^iiii 2 i # >) 



(24) 



— i= , 1 _, is the normalization fa 
by \H X ) 



where —7— , = is the normalization factor. The color part of the first trace term in the physical state is denoted 



\H 1 ) = \Tr(QG 1 G 2 )) 



where \H X ) has the order (gig 2 )- We can make the state \H X ) normalized: 



\hi) = 1 _ (26) 

y/Nf - 2N C + T/Nc 



G 



From now on, for the sake of convenience, we will simply call the color part of a trace term in the full physical state 
( f23| ) a trace state. So the projection of the state \si) on the normalized trace state \hi) is: 

(*i I hi) - ^ (27) 
We should compare the above result with that of two singlet chain states = — = |l 9l l 12 l2q) and I/2) = 



/ NS 1 1<?2 121 liq) which has the order (qgig2q) and {qg2giq) respectively. Their projections on the normalized trace 
state \h\) are given by: 

</ 2 \hi) = ^j ^ 

Hence we conclude that the CSGS accounts for about 1/N 2 of the trace state \hi) and there is a small portion (1/N^) 
of I/2) in \h\) due to the different gluon order between the two states. 

The projection of \s\) on another trace term \H2) = \Tr{QG2Gi)) is the same as on \Hi): 

(«i I Hi) = (ai I H 2 ) (29) 

since there is no order for two gluons in \H\) and \H 2 ) : \Tr{GiGi)) = \Tr{G%G\)). So the projection of \s\) on the 
full physical state \H C ) is given by: 



{SI I H c ) = ^KZU^ {D 12 + D 21 } {m) 



The probability for the parton system qqgig2 to be in \si) is obtained by: 

Prob(\ S i)) = ±-J<m\{ Sl \H c )\ 2 



±.^1 J dn(\D^\ 2 + \D^\ 2 + 2Re(D^D^)) 



(31) 



where 00 = CTtree(e + e~ — > 995132) is the tree level cross section for e + e~ — * qqg\g2- 

There is only one CSGS for qqg\g2- As the number of gluons increases, more such states arise. For the second 
example, we investigate the process: e + e~ — * qqgig2g3- The full physical color state for the process is: 

= (vfe^ 123 I^W^Ga)) + \Tr(QG 1 G 3 G 2 )) + D 213 \Tr(QG 2 G 1 G 3 )) 
+D 231 \Tr(QG 2 G 3 G 1 )) + D 312 \Tr{QG 3 G 1 G 2 )) + D 321 \Tr(QG 3 G 2 Gi))] 

Without loss of generality, we choose one trace term, e.g. \H\) = |Tr(QGiG2G 3 )) to go through our discussion. The 
other situations can be instantly obtained by permuting the gluon labels. 

As in the case of e + e _ — > qqg\g2 , when N c is large, the normalized state \hi) is approximately the same as that of 
the singlet chain state = (1/N 2 ) | lqi I12 123 l3q) , noting that has the same order of gluons as \h%). There are 
five CSGSs for the system qqg\g2g 3 '- 

{\ Si ) , i = 1,2,3,4,5} 

= {A |Tr(GiG 2 )l, 3 W , A |Tr(G 2 G 3 )l gl l 15 ) , A l^dG,)^ W , (33) 
B \Tr(GiG 2 G 3 )l^) , B |rr(G 1 G 3 G 2 )l 9? >} 

where A and B are normalization factors and given by: 

A = j , B = 1 — (34) 

N c ^Nf^l y/N* - ZN 2 + 2 

These singlets are all made up of two sub-singlets. For the first three states(i=l,2,3), one sub-singlet is formed by 
two gluons, the other is by the remaining gluon with the quark-antiquark pair, while for the last two states(i=4,5), 
all three gluons form one sub-singlet and the quark-antiquark pair forms the other one. Their projections on the 
normalized trace state \h\) are: 
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|(S1|M| 2 = I<S2|M| 2 ^W C 2 

|( S3 |M 2 ^i/^ c 6 m] 
|( S4 |MI 2 -i/^ c 2 1 j 

Ksnlhjf-l/N? 

We see that the projection squares of \si), \s 2 ) and | S4) on \h\) are all about 1/iV, 2 to the leading order in A c , while 
those of I S3) and |sg) are suppressed by an additional factor l/N^. Recalling (p3[), we may write down the parton 
orders for all these CSGSs below: 

(12)(03) -> | Sl ) 

(01) (23) - \s 2 ) 

(02) (13) - |« 3 > (36) 
(0)(123) -> | S4 ) 

(0)(132) -> |s 5 ) 

where a parenthesis denotes a sub-singlet; stands for the quark-antiquark pair and 1,2,3 for gluons. We compare 
these orders with the order in \H\): (0123). For \s\), the orders in both sub-singlets are the same as (0123), while 
for, e.g. 1^3) , the order in both sub-singlets are different from (0123) because both and 1 are not adjacent to 2 
and 3 in (0123) respectively. We conclude that the parton orders in \s\), \s 2 ) and \s4) have the least differences from 
(0123), while those in \ss) and | S5} are more different from (0123) than the former three states. In other words, each 
sub-singlet of \s±), \s 2 ) and |s<j) conserves the parton order in the trace state at most, so they have the largest overlap 
with the trace state. The reason that |s 3 ) and |sg) are more suppressed is that their sub-singlets do not all keep the 
parton order in the trace state. The inner product of any two states in |si), \s 2 ) and \s&) is about l/N%, while that 
of one state in |si), \s 2 ) and \s^) with one of \ss} and |sg) is about or higher. 

This rule can be further verified for e + e~ — + <?95i92.9394- CSGSs with the leading contribution to \Hi) = 
\Tr(QGiG 2 G 3 G±)) are: 

{(123) (04), (01)(234), (12)(034), (23) (014), (012) (34), (0)(1234)} (37) 

Their projection squares on the normalized trace state \hi) are all 1/-/V 2 to the leading order in 1/N C . The inner 
product of any two different states in the above list is about 1/N%. Note that we only include states with two sub- 
singlets because the projection squares on \h\) for those with three or more sub-singlets are at least suppressed by 
additional 1/A C 2 as compared to those with only two. 

For cases of n-gluons where n < 6, the above conclusions can be verified by direct calculation. For cases of n > 6, a 
formal proof is given below. Suppose any two different leading CSGSs have an inner product of 1/N% to the leading 
order in \/N c in n-gluons case. Let us prove that the same statement holds for (n + l)-gluons case. We focus on two 
different CSGSs which we call A and B: 

A : (/)(//) = (0, 1, 2, • • • , ki, fci + h + 1, • • • ,n + l)(fc x + 1, ki + 2, • • • , ki + h) 

B : (IH)(IV) = (0, 1, 2, • • • , k 2 , k 2 + l 2 + 1, • • • , n + l)(k 2 + 1, k 2 + 2, • • • , k 2 + l 2 ) [M) 

where stands for qq and the parentheses for sub-singlets. Each CSGS has a gluonic sub-singlet with li and l 2 gluons, 
respectively. For n > 6, at least a common piece of [01,02,03] in two CSGSs must exist. Here a\, a 2 and 03 denote 
three consecutive gluons in a sub-singlet of one CSGS. When we calculate the inner product of A and B, we get a 
factor of N% to the leading order in N c from this common piece. Corresponding to A and B in (n + l)-gluons case, 
there exist two CSGSs A' and B' in n-gluons case where all parts of A'/B' are the same as A/B except that the 
common piece [01,02,03] is changed to [01,02] (03 as the new added gluon to A and B) or [02,03] (ai as the new 
added gluon to A and B) or [01,03] (02 as the new added gluon to A and B). Note that any one of three cases are 
equivalent to our proof. The contribution from the common piece, say [ai, 02], to the inner product of A' and B' is N c 
to the leading order in 1/N C , while that from the common piece [ai, 02, 03] to the inner product of A and B is about 
N%. The normalization factor of about N c (™ +2 ) for (n + l)-gluons case and Nc^^ 1 ^ for n-gluons one must also be 
taken into account. Considering contributions from both the common piece and the normalization factor, we find that 
the inner product of A and B is the same as that of A' and B' to the leading order in l/N c . The same strategy can be 
also used to prove that the inner product between any leading CSGS and its corresponding trace state is about l/N c . 
Hence we summarize the following rules for the CSGS with the leading contribution to a trace state: 

• The state has only two sub-singlets to the leading order in 1/N C . Its projection square on the trace state is 
about l/iV c 2 . 

• Each sub-singlet keeps its parton order as that in the trace state. 
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• The inner product of any two different leading contribution states is about 1/N%, while that of a leading state 
with a non-leading one is about 1/iVj or higher. 

According to the above rules, we can determine for the process e + e~ — > qqgig2 • ■ ■ Qn-\9n the total number of 
CSGSs with the leading contribution to a trace state. We suppose this trace state corresponds to the trace term 
\Tr{QG\G2Gz ■ ■ ■ G n -xG n )). Hence, for the configuration (qqgig2 • ■ • gn-i9n)i there are (n — 1) states with one sub- 
singlet formed by 2 gluons and the other by the remaining gluons with the quark- antiquark pair. Obviously, there are 
(n — k + 1) states with one fc-gluon sub-singlet and the other formed by the remaining partons. Note that all these 
sub-singlets keep their parton orders as that in the trace state. Altogether, the total number is: 

1 + 2 + 3H hn- 1 = -n(n- 1) (39) 

We hereafter denote this number by m = \n{n— 1). 

Now we start to discuss how to estimate the total probability of the CSGS. One problem we are facing is that CSGSs 
are not orthogonal to each other. One cannot obtain the total probability by summing all individual probability values 
for each state. The non-orthogonality would make the total probability derived in this way different from its correct 
value because the contribution from the overlap of any two states would be counted multiply. So for non-orthogonal 
states, the orthogonalization is needed. Let us illustrate the procedure more explicitly. If we have a set of non- 
orthogonal states: 

{\ f i),i= 1,2,- --1} (40) 
suppose we have finished the orthogonalization and find a set of orthogonal states as follows: 

{]/;>,* = i,2,- ••/} (4i) 

The total probability for the set of states is: 



Prob(J2\fl)) = i-EiIdn\(fi\H c 



<?0 

£ if dn IEp ° P (fl I Tr{QG P{x) G Pm ■ ■ ■ G P(n) )\ 



(42) 



where oo = Ct ree (e + e~ — * qqgig2 ■ ■ ■ <7n) is the tree level cross section for e + e~ — > qqgig2 ■ ■ ■ 

As stated above, it is necessary to orthogonalize the non-orthogonal CSGSs to give their total probability. But 
actually this job proves to be extremely difficult. The greatest obstacle lies in that the number of CSGSs (which is also 
the dimension of the transformation matrix from original states to orthogonal ones) is huge when the gluon number is 
large, and finding the transformation matrix with this huge dimension is almost impossible. All possible CSGSs must 
be taken into account, not only those with the leading contribution to a specific trace term of the full physical state 
\H C ). The reason for this is simple: \H C ) includes all trace terms each of which corresponds to one order of gluons, 
hence a non-leading singlet for one trace term is not necessarily still the non-leading one for other trace terms with 
different gluon orders. Another big obstacle is that calculating the momentum function D p is a very complicated 
task, especially for multigluon processes. To some degree, these two difficulties make the orthogonalization of more 
theoretical value. 

The orthogonalization is, however, really simple and heuristic in some cases. Let us look at such a case. Assume 
that the physical state is fixed at one trace term, e.g. \H) = \Tr(QG\G2Gz ■ ■ ■ G n -iG n )), we now try to estimate the 
probability of the CSGS in this case. Before we orthogonalize these singlets, we make some further approximations: 

• We only take into account the leading contribution singlets as determined by the rule on page ||. All the inner 
products of two different leading singlets are supposed to be the same l/N^. In fact, they are the same to the 
leading order in 1/N C . 

• There is freedom to choose the transformation matrix (the definition is given below) . Considering the symmetry 
of inner products assumed in the previous item, we choose a simplest form for the transformation matrix: all 
the diagonal elements are assumed to be equal, and so are the non-diagonal elements. 

Keeping these approximations in mind, we will proceed with the problem. Now we have m non-orthogonal states: 

{|fli),i=l,2,-..,m} (43) 
where the inner product of any two states is given by: 
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Our goal is to find to orthogonal states based on (|4^). These states are denoted by: 

{|s<> , » = l,2,...,m} 
They are related to original states ( |43| ) by a linear transformation: 

\s'i) = Ua \8j) (45) 

where the transformation matrix U is a real to x to matrix. It is not the only transformation which transforms the 
original non-orthogonal states to orthogonal ones. Different transformation matrices correspond to different orthogonal 
sets. Two transformation matrices are connected through a unitary transformation, but the transformation matrix 
itself is not unitary because it associates orthogonal states with non-orthogonal ones. Any transformation matrix is 
equivalent to others in calculating the probability. The inner product of two new states is: 

(4 | 4) = S(i,j) = UikUjiisk \si}=J2 UikUjiCki (46) 
k,i k,i 

As we state in the above, from the symmetric property of Cki in (0), we choose a symmetric matrix U . We search 
for a solution for which all diagonal elements have the same value and also all non-diagonal elements are equal to 
each other: 



Ua = a 

Uij = b, for i ^ j 



where a and b are real numbers. Thus (^6|) becomes: 



J2k,i UikUuCki — i 

J2k i UikUjiCki = , for i ^ j 



Substituting @ into @, we obtain: 



a 2 + (m - l)b 2 + jfs(m - l)ab+ ^{m - l)(m- 2)6 2 = 1 

2ab + (m - 2)b 2 + ^a 2 + ^(m - 2)ab +^{m- 2)(m - 3)b 2 = 



(47) 



(48) 



(49) 



We find the exact solution for (|49|): 

N 2 



Q {m 2 -'im + 2 + (m-l)N 2 )y 2 + 2{m-l)y + N 2 /' h W V ' A)) 



where y is given by: 



_ -(m + N 2 -2)± y/(N 2 + 1)to + N± - 2N 2 -2 
V (to ,- 2)(m + N 2 - 3) 1 ' 

We note that this solution is finite when to — > 2, and therefore regular for all positive values of m. We choose the 
solution with a positive sign, as this corresponds to a unit matrix, U = 1, in the limit N c — > oo. The orthogonal 
states are: 

\s' i ) = a\s i )+bJ2\s j ) (52) 
The overlap between the trace state \h) and the states \sj) is given by 

Ei< s ^)i 2 = ^( a + ( w - 1 ) 6 ) 2 ( 53 ) 

i c 

When N c is large we find y ~ ~2W* anc ^ ° 2 — 1- Thus to order 1/N 2 the probability behaves as which increases 
monotonously with increasing number of giuons. To estimate the probability we introduce the singlet chain state |/) 
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into the list of states to be orthogonalised. Note that |/) has the same gluon order as the chosen trace state \h) to 
which the CSGS {|sj)} or {|s^)} correspond. We define the new orthogonal states as follows: 

|<}=aK)+/3|/) (u] 
\f')=l\f)+SJ2 k \s' k ) 

When N c is large, keeping only the leading terms in N c , (H^) becomes 

| s '.)~| Si )--L_£| Si ) (55) 

C 3& 



which implies (/|s£) ~ l/N c . a, /3, 7, <5 are four parameters to be determined from the following orthogonality: 

(56) 



( s "\ s i) = 1 

(s'l\s'!) = for i± j 
(/'I/') = 1 

(/X) = o 



From the first and second equation of (p6|), we obtain a = 1 and (3 = —2/N c . From the last two equations of (p6|), we 
get to order 1/N% 

In the large N c limit, the squared projections of |/') and \s") on the corresponding normalized trace state \h), i.e. the 
probabilities of |/') and \s") in \h), are given by: 



Pf>*\(f'\h)\ = 7 2 (l + m/^ c 2 ) 2 -i^ 



(58) 



In the last aproximations we have added terms of higher order in 1 /N£ in such a way that the total probability adds 
up to one for all values of m. Thus the total probability for all \s") states becomes mP s » ~ i+j^/N' 1 ' which increases 
monotonously with growing gluon number (or to). Correspondingly Pfi decreases monotonously with increasing m. 

Summarily, what we have derived gives us some valuable messages about the properties of the CSGS. Here are some 
of them: (a) The probability of color separate singlets cannot be neglected. They may account for a considerable 
fraction of the total events. We try to estimate the probability to produce a CSGS, disregarding any dynamic effects 
which might disfavor its hadronization as a separate unit, (b) The set of states corresponding to different string 
configurations must necessarily be orthogonal, as they correspond to different values for observables. In the large N c 
limit, CSGSs do not appear. A trace state is just the singlet chain state having the same gluon order. All singlet chain 
states are orthogonal to CSGSs. For finite N c , neither are the three sets of states, trace states, singlet chain states 
and color separate states, orthogonal among themselves, nor is a set of states orthogonal to any other. These states 
are not possible candidates to represent pure string states. Keeping only the first correction term in an expansion 
in 1/N%, it is, however, possible to construct modified orthogonal states which corresponds to a specific trace state 
\hi): {/■, {sy}}- Thus they form an orthogonal set of states, which may be associated with specified string (or chain) 
configurations. Here we propose a first estimate for the group weight corresponding to a specific string configuration 
obtained by projecting these string candidates upon a specific \hi). (c) For a trace state, leading CSGSs obey the 
rules on page ^|. In particular the rules are crucial for us to build a practical model to describe the production of color 
separate states. This is the subject which we will discuss in the next section. 



V. A COLOR SEPARATE STATE MODEL 



According to the arguments in the last section, the total percentage of color separate states should not be negligibly 
small, especially in multigluon process. It is of significance to introduce CS states in the current event generators to 
see how they influence the event shape and properties. In this section, we will build such a practical model based on 
and implementing the rules which we derived in the former section. 

Firstly, the difference between two terms, the CSGS and CS state, should be clarified. We mention in the in- 
troduction that to describe the interface between the perturbative stable partonic state and the hadronic one, it is 
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not enough to know the momenta and colors of all partons, we also have to know how they are connected. Their 
connection can be described by a string or cluster chain configuration where each color charge must be associated 
with a partner anticharge. The connection or the string/chain ordering cannot be specified by perturbative QCD 
alone. It depends on the nonpcrturbative, soft, properties of the strong interaction. If the string or cluster chain state 
contains a closed loop formed by gluons we call it a CS gluonic sub-singlet which is a non-perturbative object, while 
a CSGS is defined by the colors of partons which is of perturbative nature. A CSGS is not necessarily a CS state. 
In this section, we use JETSET to implement the fragmentation which is a non-perturbative soft process. Hence we 
start using the term "CS state" in most places to describe the true objects formed at the end of the parton cascade 
or at the beginning of the hadronization. 
The outline of the model is as follows. 

1. We use a parameter R to describe the total percentage of CS states in all events. We use JETSET to produce 
a configuration of partons where for the singlet or for each sub-singlet of the type qg\g2 ■ ■ ■ giq, the color flow 
starts at the q, connecting gluons one by one and ends at q. This color connection corresponds to a trace term 
in the physical color state \H C ) which we call a trace state. 

2. We only consider those CSGSs with leading contributions to the trace state. For a configuration qgig2 • ■ ■ giq, 
there are m = |Z(i — 1) leading states according to (|39|), Among them there are (I — 1) CSGSs of 2-gluon type, 
(1—2) states of 3-gluon type,...,(i— fc+1) states of fc-gluon type, etc.. Any two different CSGSs have approximately 
the same inner product(~ 1/N%). Each state's square of projection on the trace state(corresponding to the 
current color configuration) is the same (~ I/AT?). So all CSGSs play an equal role. However, not all of these 
states finally contribute to the event. Thanks to the fact that JETSET records the whole branching history 
for all partons, we can exclude those CSGSs where all gluons in the gluon-type sub-singlet are emitted from a 
single octet gluon. After filtered by this procedure, suppose there are m' leading CSGSs remaining. 

3. We consider two kinds of weights to select a CS state corresponding to the given configuration. One is the 
constant weight, i.e. each CS state is selected with equal probability. This is a natural choice if we only consider 
the color part and neglect the momentum or phase space sector. The other is T(or A) measure weight(we 
will come to this topic later in details). For each CS state, we have a T- weight. So there are m' T- weights 
corresponding to m' CS states. A CS state for the current configuration is selected according to these weights. 

4. The fragmentation of the selected CS state is implemented by the Lund string model. The difference between 
an ordinary state of the type qg\gi • • ■ giq and a CS state is that the former state forms an open string while the 
latter forms a closed string and an open one. 

In order to better understand our Monte Carlo simulation results, we should use a little more space for illustrating 
the physical meaning of the T-measure. We will address two subjects: (a) How is the T-measure defined and what's 
its physical meaning? (b) Its relationship to the A-measure in the Dipole Cascade Model(DCM). 

For a qq system the produced hadrons are evenly distributed in rapidity, which means that their energy-momentum 
four-vectors, when plotted one by one, are distributed around a hyperbola as seen in Fig.|l|a. For a qqg system, 
the two string pieces will produce hadrons such that the momentum four vectors lie around two hyperbolae(Fig.|l|b). 
This corresponds to three jets along the parton momenta, and the jets are smoothly connected to each other. For a 
multigluon system it is possible to generalize the hyperbolae in the qq and qqg cases and define a timelike curve in 
energy-momentum space. This curve (which we call the x-curve) follows the color ordered parton momenta in such a 
way that the corners are smoothed out with a resolution power given by a parameter mo- The directrix curve for 
the string parton state is obtained by placing the four-momenta of all the partons one by one according to their color 
ordering, see Fig.|l|c. When mo — > 0, the x-curve becomes the directrix. Let £ be a parameter which runs along the 
directrix. To every point on the directrix there is a corresponding point on the x-curve. So the x-curve which 
is determined through the differential equations |llj: 

dx,(0 = ^(q-dA)q, 

x^0) = Ajp) = 1 ' 

where q^) = — 2^(0 is tangent to the x-curve. T is defined by: 

d(lnT) = ^ = #2 = dX 



o y/q" (60) 

T = exp(X) 

Note that T is just the exponential of the area spanned by the x-curve and the directrix. T can also be interpreted as 
the exponential of the invariant length of the x-curve which we call the A-measure. The square amplitude to obtain 
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a certain string state with n partons is given by T ^ k\ ■■■ k ) wnere ki ? k 2 , • • • , k n are four- momenta of n partons. 
There is a recursive expression for T n : 

= 7i+i9j + |(1 + 7i+i)*;+i 
7i+ i = l/[l + ^i] (61) 

T^ 1 = 7273 ■■■In 

The T-measure defined as above depends on the initial value of q^. It is possible to start with the initial value qo = 
for an open string. We regard the sequence fci, k 2 , ■ ■ ■ , k n as the first half of a full directrix period, the second half 
is k n , k n -i, • ■ ■ ,ki, which is in the opposite order. This is one period: fci, fc 2 , • • • , k n , k n , k n -i, • ■ • , k\. After several 
periods of recursive calculation according to ( fill) > we will get for a period a stable T which is insensitive to the initial 
value go • For a closed string where there is no obvious beginning, as a choice, we also use this strategy to get a stable 
T. 

In the string fragmentation model, the A-measure is the effective rapidity range given by: A ~ ln[(ki + ki + i) 2 /m^] 
Hence, according to the weight defined by ^ = exp(— A), we may tell which string parton state is "short" and which 
is "long". Here a "short" string means the string with small A. The weight implies that a short string parton state is 
more favored than a long one. 

For a configuration (qgig2 ■ • ■ dnQ) produced by JETSET, suppose an allowed CS state is made up of a separate 
sub-singlet (<?i<?i+i ■ ■ ■ gi+j-i) of the j-gluon type and a sub-singlet [qg\ ■ ■ ■ gi-igt+jgi+j+i • ■ ■ g n q) of the ordinary type. 
The total weight is: 

i_ 1 

T T 1 (k i ,ki + i,---,k i+ j- 1 ,k i )T2{k q ,k 1 ,---,ki- 1 ,k i+ j,k i+ j +1 ,--,k n ,k^) (62) 

= exp(—\i - A 2 ) 

There is a difference between T\ and T 2 : T\ is the T-measure for a closed string, so it is written as: 

Ti = 7fc I+1 7fc l+2 • ■■lk i+j -i7k i = exp(-Xi) (63) 

where the recursive sequence (fci, fci+i, • • • , fci+j_i, k{) starts at fci to fci + j_i and then back to fci; and Ai is approxi- 
mately: 

Ai ~ ln(s iti+ i/ml) + ln(s i+ i 4+2 /ml) H h ln(s i+ j- 2 ,i+j-i/ml) + ln(s i+J ^i :i /ml) (64) 

where Sij = (fci + kj) 2 and the last term is the additional term compared to an open string. There is no such a term 
for an open string formed by the same sequence of momenta. We also use the A-measure for the full configuration 
(?5i32 • • • ffn?) without interruption by separate sub-singlets: 

A ~ ln(s lt2 /ml) + ln(s 2 ,3/ml) H h in(s n _ 2 ,n-i/ m o) + M s n-W TO o) ( 65 ) 

We can compare Ao with Ai + A2 by investigating their difference (Ai + A 2 ) — Aq: 

Ai + A 2 — Ao 

~ ln(s i+ j-^i) + ln(si^i ti+J ) - ln(si-i t i) - ln(s i+J ^i^ i+j ) (55^ 

_ , / (fc i+J -_i-fc i )(fc i _i-fc i+j ) s 

- ir ^(k i - 1 -k i Kk i+j - 1 -k i+j )J 

where fci_i and fci, fci +J _i and ki + j are momenta of two gluons in neighborhood in the color flow respectively, while 
fci + j_i and fci, fci_i and fci+j are not adjacent. Normally, we have: 

(ki+j-i ■ ki){h-i ■ k i+ j) > (ki-i ■ ki)(ki+j-i ■ h+j) (67) 

which implies Ai + A 2 > Aq. That is to say: the length of a CS state is usually longer than the ordinary state without 
interruption. We expect that the shorter the length of a string state is, the more favored it would be. We can find 
some configurations for a CS state where its length Ai + A 2 is shorter. It is conceivable that these configurations are 
more favored and have larger production weights. For example, if all the j gluons in the separate sub-singlet plus two 
neighbor gluons to its two ends are soft, then the weight must be large because separating the j-gluon as a sub-singlet 
will not change the original Ao very much. 

To see the situation more easily and intuitively, we resort to the A diagram which is often used in the DCM. Fig.|| 
is the diagram of the effective rapidity range (i.e. A) for a parton configuration pQ ]. In this configuration, there 
are nine gluons. The order along the color flow is shown in the figure which reads qg\g 2 ■ ■ ■ 399- The length of the 
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folded line ql23456789<7 is Ao- The length of the folded line between two neighboring gluons, e.g. i and i + 1, is 
Zn(s, l i+i/mo)- In the figure, there are 11 ends labeled by q, q and nine gluons respectively. We call these ends twigs. 
Along the line qq, there are 5 branches. In the branch which starts at D, there are 4 twigs, two of which, 6 and 7, 
are from the sub-branch starting at F. For this parton configuration, we will look at two examples of CS states to 
see how their A measures differ from Ao- The first one is: ((?123g)(456789), which is shown in Fig.pk. The length of 
the dashed line shows the A surplus compared to Ao- The surplus is denoted by <5A: 8X = 2\DH\. The second CS 
state is: (gl234569g)(78) shown in Fig.|b. The A surplus is: SX = 2(\FG\ + \DF\ + \HI\ + \DH\). Of course, the 
first CS state is more favored than the second one though it has more gluons in its gluon sub-singlet. If we look at 
them more carefully, we will find that the gluon sub-singlet of the first CS state is made up of all gluons from two 
adjacent branches, while the sub-singlet of the second is formed by part of gluons in the two branches. We also see 
that the shorter the line \DH\ is, the larger production weights the CS states have. Here \DH\ is just the rapidity 
difference between the two adjacent branches. We see that those CS states with their gluon sub-singlets formed by 
adjacent twigs (note that twig is the smallest branch) have large production rate. In this case most gluons in gluon 
sub-singlets are soft ones. 

We summarize the above arguments that most of CS states produced in the process have the following configurations: 

• The gluon sub-singlet of a CS state is made up of a series of neighboring gluons. If part of a branch joins the 
sub-singlet, the rest of this branch should also join it. Any separation of one part from the rest of the same 
branch would result in the suppression of the production weight. The whole branch is like a skeleton. The 
principle says, if the skeleton is in, the whole body should be in. 

• The rapidity difference between two branches at the two ends of the continuous branch sequence is as small as 
possible. This implies that all branches are as far from the root (or beginning) of the cascade tree as possible. 
Hence those CS states with their gluon sub-singlets formed by adjacent twigs have large production rate. 

Now we begin to discuss our Monte Carlo results which are shown in Fig.^-^. Fig.|| and || show how CS states 
influence the global properties of events, while Fig.|](a-d) and 0(a-d) shows other properties in the case of the T- weight 
and the constant weight respectively: the distribution of the gluon number in the gluonic sub-singlet of a CS state, 
the multiplicity distribution in the whole and the restricted phase space, and the rapidity distribution for c-events in 
the central region. The constant weight is defined as that each possible CS state is chosen with the same probability, 
and the T-weight is discussed earlier in this section. We show three groups of results: no CS states, 100% and 30% 
CS states. These results show that for the situation of the constant weight, there is observable though small deviation 
between the result with CS states and that without CS states, while for the situation of the T-weight, the deviation is 
even smaller. One may easily accredit the small deviation in event properties to the effect of the CS state. This 

is not truely the case. Here two questions should be distinguished: (a) How large is the pobability for CS states to 
appear? and (b) How large is the deviation between the hadronic events with and without CS states included? They 
are two distinct concepts by nature. We recall that there are many non-orthogonal leading CS states corresponding 
to an ordered trace state. Though a single leading CS state contribute only l/N^, the total contribution may actually 
be large for many CS states. In fact we cannot determine the total probability at the perturbation level as shown 
earlier in the former section. Hence we include a parameter R to denote the total probability of CS states which 
we will give results in two cases R=30Note that all the above statements are refer to question (a). But when we 
talk about the small deviation between events with and without CS states, we are addressing question (b). One can 
understand this deviation as the result of the soft gluon dominance in the cascade and the infared property of global 
observables, considering that the global shapes are mainly determined by the cascade instead of the hadronization 
and that the effect caused by isolating a group of gluons (mostly soft continuous gluons in the color sequence) should 
not be large. We know earlier in this section and later in the following paragraph that T-weight favors the CS states 
with softer and less gluons in their gluonic subsinglets which bring little change to global and local observables, while 
the constant weight gives each CS state equal probability which leads to larger contribution from those CS states with 
harder and more gluons in their gluonic subsinglets whose observables may differ significantly from those conventional 
events without CS states. So on average the constant weight gives larger difference from conventional events than the 
T-weight. 

It is meaningful to investigate local observables especially in carefully chosen phase space regions and types of 
events, where the difference between CS events and Non-CS ones may be more clearly seen. Fig.^a shows the 
probability distribution of gluon number in the gluon type sub-singlet of a CS state. The T-weight sampling is 
used to select a CS state for a given parton configuration. The distribution shows a perfect exponential decreasing 
feature: P{n g ) ~ exp(— 0.9n g ). The mean value of the gluon number < n g >^ 2.9. Fig.^a is the n g distribution 
with the constant weight. Its decreasing behavior is much slower than that with the T-weight. It has the mean value 
< n g >~ 3.5. We know that the T-weight favors the CS state with soft gluons in its gluonic subsinglet. More gluons 
are in its gluonic subsinglet, more probably it gives a smaller T-weight. Hence the T-weight favors the smaller number 
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of gluons in the gluonic subsinglet of a CS state. Therefore the < n g > distribution with the T-weight is steeper than 
that with the constant weight. FigJ^D is the multiplicity distribution for charged particles using the T-weight. The 
solid line is the result without CS states, and the dashed line is that CS states are produced with the total fraction of 
30%. We see that there is little difference after introducing CS states. It is not surprising because we know that in most 
cases the gluonic sub-singlet of a CS state is composed of soft gluons, which brings little change to the multiplicity of 
the original string connection. Fig.^b is the multiplicity distribution for charged particles in the constant weight. An 
apparent deviation from that without CS states can be seen. The broader distribution gives a larger mean multiplicity 
than in the T-weight case. This phenomenon is consistent with our above analysis that compared to the T-weight 
case, the constant weight enhances the probability of CS states with harder gluonic subsinglets which give rise to 
larger multiplicity. Fig|j|c is the multiplicity distribution of charged particles in the T-weight in the rapidity range 
— 1 < y < for the c-events. The rapidity axis is chosen along the added cc momentum. The selection criterion is 
that the angle of cc is smaller than 110°. The solid line is the result without CS states, and the dashed and dotted 
line are that CS states are produced with total fractions of 100% and 30%, respectively. The reason for choosing the 
c-event is that for these events it is easy to tag the primary quarks. Normally there should be a rapidity gap caused 
by CS states between the cc cone and the cone in the opposite direction. That is to say, in the rapidity region around 
y ~ 0, there are fewer particles than the case without CS states. This effect can be seen in the figure. The probability 
to have fewer than 2 charged particles in the range — 1 < y < in the case with CS states produced is obviously larger 
than that without them. This is just what we expect a sign (though a weak one) of a depleting area or a rapidity 
gap. Fig.fjjc is the result with the constant weight. No sign (even a weak one) of a depleting area is found. This is 
possibly due to the broader multiplicity distribution brought by the constant weighting (see also Figj^b). Fig.||d is 
the rapidity distribution of charged particles for the c-event in the T-weight. The rapidity axis is chosen along the 
added cc momenta. The selection criterion is that the angle of cc is smaller than 110°. The solid line is the result 
without CS states, and the dashed line is that CS states are produced with total fraction of 30%. The difference of 
the largest and the most interest is in the region —2 < y < 0, where the distribution with CS states is a little lower 
than that without them. This is a weak indication for the particle defeciency in the central region. This phenomenon 
is absent for Fig.[j]d where the constant weight is used. On the contrary, the central region of the distribution with 
CS states is clearly higher than that without them, which implies a more populated central area exists for the case of 
CS states. 

In the literature, there are a number of reconnection models which can describe the production of CS states in 
a variety of reactions [f7||8||l^]. Ref. Jl2| gives a reconnection model(which we call RMDCM for short) which is 
incorporated into the dipole cascade Monte Carlo ARIADNE. In some sense, its strategy of color reconnection is 
similar to ours, but it is more complicated and more phenomenologically based. Our model provide a significantly 
simpler recipe for reconnection. The results for e + e~ annihilation produced by our model and RMDCM are similar. 

Noticing that the CS states with soft-gluon sub-singlet have large T-weights and hence have large production rates, 
the results of our model have some similarities with those from the Soft Color Interaction (SCI) model pi] . We will 
in the future study implications from our model for deep inelastic scattering. 

VI. SUMMARY 

We use the method of color effective Hamiltonian H c to study properties of color separate states in e + e~ — » qq + ng. 
The full physical color state is \H C ). There are n! trace terms in \H C ), each corresponds to an order of n gluons. We 
call each trace term a trace state. Up to 0(1/ N^), there are m = ^n(n — 1) CSGSs with the leading contribution to 
a trace state. We obtain the following rules for a leading CSGS: Up to the leading order in 1/N C , (a) The state has 
only two sub-singlets. Each sub-singlet keeps its parton order as that in the trace state at most, (b) Its projection 
square on the trace state is 1/N^. (c) The inner product of any two different leading states is 1/N%, while that of a 
leading state with a non-leading one is 1/iV? or higher. 

Due to the non-orthogonality of CSGSs, it is necessary to orthogonalize them to give the total probability. Actually 
the general orthogonalization proves to be extremely difficult. But it is simple and heuristic in some circumstances, 
e.g. the physical state is fixed at one trace state. In this case, we can estimate the probability of the CSGS 
through orthogonalization under some further approximations: (a) We take into account leading CSGSs and the 
singlet chain state with the same gluon order as that of the trace state, (b) We choose a simple form for the 
transformation matrix. The orthogonalization procedure has two steps. The first step is to find the orthogonal states 
{|sj),i = l,---,m} from m leading CSGSs {|sj),i = 1, •••,m}. The second step is to introduce the singlet chain 
state {|/)} which is approximately the trace state {\h)} into the list, and to obtain the final set of orthogonal states: 
{|/'), {|s"), i = 1, • ■ • , m}}. In the large N c limit, the probabilities of \ f) and one \s") in the corresponding normalized 
trace state \h) are given by: Pf ~ i+m/iv 2 and P s // ~ (1 — 1+r ^/ N i )/m — 1/N%. The total probability for all \s") 
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states are mP s " ~ 1 — 1+r ^/ N -2 — m/N%. We see that Pfi decreases while mP s " increases monotonously with the 
growing gluon number(or m), and Pfi +mP s » conserves the total probability 1. These results are for large N c . When 
N c is finite, we always span the pole at m = N%/2 for P s " which is unphysical. Hence the present result is invalid for 
finite N c . In this case, We have to consider all non leading terms in N c which is too complicated to have an explicit 
result. 

The rules for the leading CSGS are crucial for us to build a practical model to describe the production of color 
separate states in event generators. Here are the outline of the model: 

1. We use a parameter R to describe the total fraction of CS states. We use JETSET to produce a parton 
configuration with a specific parton connection. We let it correspond to a trace state. 

2. We only consider CSGSs with the leading contribution to the trace state. According to the cascade history 
recorded by JETSET for the current configuration, we can further exclude CS states whose gluon-type sub- 
singlets are emitted from a single octet gluon. 

3. We consider two kinds of weights in selecting individual CS state for an event. One is the constant weight. This 
is a natural choice if we only consider the color part and neglect the momentum or phase space sector. The 
other is the T(or A) measure weight. 

4. The fragmentation of the selected CS state is implemented by the Lund string model. Note that the gluon 
sub-singlet forms a closed string whose fragmentation should be managed with care. 

If we use T measure weight, we find that most of CS states produced in the process have the following configurations: 
(a) The gluon-type sub-singlet of a CS state is made up of all gluons in several continuous branches. If part of a 
branch joins the sub-singlet, other part of this branch should also join it. Any separation of one part from the rest 
of the same branch would result in the suppression of the production weight, (b) CS states with their gluon-type 
sub-singlets formed by adjacent twigs of the branching tree have large production rate. From these points, we see 
that, CS states whose gluon-type sub-singlets consist of soft adjacent gluons have relatively larger T-weights than 
other types of CS states and therefore have larger production rates. 

The comparisons of global and local event properties between events with and without CS states in the T- and the 
constant weight are given. Generally the difference of results between two cases, those with and without CS states, 
is found to be small, though we find the deviation in the constant weight is a little larger than that in the T- weight, 
and in the central phase space area a weak sign of the rapidity gap is shown in the T-weight case. The main reason 
lies in the fact that soft gluons dominate in the cascade. Considering the infared safe properties for all observables, 
the effect caused by isolating a group of gluons (mostly soft continuous gluons in the color sequence) should not be 
large. T-weight favors CS states with soft gluonic sub-singlets which should not bring significant change to global and 
local observables, while the constant weight gives each CS state equal probability which enhances contributions from 
those CS states with harder gluons in their sub-singlets whose observables may differ significantly from conventional 
states without separate subsinglets. So on average the constant weight gives rise to larger deviation from conventional 
events than the T-weight. Conclusively we find that in the high energy electron-positron annihilation, global and 
local properties are mainly determined by the momentum configuration of final partons, and are almost insensitive 
to the way of color connection. 
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(a) (b) (c) 

FIG. 1. (a) In a qq system, the hadron momenta are distributed around a hyperbola in energy- momentum space, (b) For a 
qqg system, the hadron momenta are distributed around two hyperbola, (c) For a multigluon state, the hadron momenta are 
distributed around a curve in energy-momentum space(called the x-curve) which smoothly follows the (color-ordered) parton 
momenta. 
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FIG. 2. Diagram of the effective rapidity range (i.e. A) for a parton configuration. 
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FIG. 3. Diagram of the effective rapidity range (i.e. A) for CS states, (a) (gl234569g)(78) (b) (ql2Sq) (456789) 
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FIG. 4. This figure is for comparing global quantities of events (sphericity, aplanarity, thrust and oblateness) among cases 
where CS states are produced with 0%, 100% and 30% probability. The CS states are chosen with constant (or equal) weight. 
The energy is 91 GeV. Solid line: no CS states; Dashed line: 100% CS states. Dotted line: 30% CS states. 
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FIG. 5. Same as Fig.H except that the CS states are chosen with T-weight. 
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FIG. 6. The energy is 91 GeV. T-weight is used to select CS states, (a) The probability distribution of gluon number in the 
gluon type sub-singlet of the CS state, (b) Multiplicity distribution of charged particles. Solid line: no CS states; Dashed line: 
30% CS states, (c) Multiplicity distribution of charged particles in the rapidity range — 1 < y < for c-events (i.e. the initial 
quark-aniquark pair produced at the electroweak vertex is cc). The rapidity axis is chosen along the added cc momentum. The 
selection criterion is that the angle of cc is smaller than 110°. Solid line: no CS states; Dashed and dotted line: 100% and 30% 
CS states, respectively, (d) Rapidity distribution of charged particles for c-events. The rapidity axis is chosen along the added 
cc momenta. The selection criterion is that the angle of cc is smaller than 110°. Solid line: no CS states; Dashed line: 30% CS 
states. 
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FIG. 7. The same as FigJq except that the constant weight is used in selecting CS states. 
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